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INTRODUCTION 


Numerical  mettiode  for  obtaining  bound-state  eolutlons  of  the  one- 
dlmenelonal  Schroedlnger  equation  can  be  broadly  classified  Into  two  types; 
numerical  Integration^*’^  (shooting  methods)  and  matrix  meAods.  The 

techniques  to  be  described  here  belong  to  the  integration  category.  They 
have  been  testod  on  a variety  of  problems  and  have  proven  to  be  efficient 
and  trouble  free  In  operation.  They  have  several  noteworthy  features. 

They  are  very  efficient  of  computer  time  (typical  run  times  will  be  given 
In  Section  IV).  No  difficulties  are  encountered  In  converging  to  the  desired 
solution  in  problems  with  double  minimum  potentials,  even  wlfo  a very  hl|h 
barrier  between  the  minima.  No  overflow  problems  are  encountered  In  the 
classically  forbidden  regions,  so  no  special  programming  precautions  have 
to  be  taken  In  this  regard.  Finally,  the  methods  to  be  presented  here  can 
be  easily  and  directly  gmerallsed  to. a practical  method  for  calculating  dm 
eigenvalues  and  eigenfunctions  of  the  multichannel  Schroedlnger  equation. 

The  calculation  of  eigenvalues  will  be  discussed  In  Section  n.  The 
numerical  techniques  to  be  Introduced  are  the  log  derivative  and  renonnal* 
laed  Numerov  meti&ods.  The  log  derivative  formulation  has  been  discussed 
previously  as  a method  for  computing  the  S-matrlx  for  scattering  prob- 
lems; foe  renormallsed  Numerov  method  Is  new.  The  calculation  of 
eigenfunctions  by  the  renormallsed  Numerov  method  will  be  discussed  in 
Section  m.  Also,  a very  useful  Interpolation  formula  for  calculating 
values  of  the  eigenfunctions  not  located  at  the  grid  points  will  be  derived. 
Finally,  foe  results  of  example  calculations  and  a general  discussion  will 
be  presented  la  Sectloa  IV. 
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more  detail  et  the  end  of  the  eectlon.  An  ennnftlon  to  the  n^thod  deecrlbeda 
above  o^cara  la  the  case  of  e-wave  bound  etatee  In  an  attractive  potential 
well  which  haa  no  abort  range  repulaioo*  caae  there  la  no  claaal> 

caUy  forbidden  region  to  the  left  of  the  potential  well  and  ao  we  nmat  aat^ 
s 0 and  uae  the  exact  boundary  condltlona. 

The  Inward  Integration  la  done  flrat  and  atopped  at  the  grid  poln^ 
where  the  relation 

17  _ 

flrat  occura.  Thla  la  our  choice  for  the  matching  point  Xj^.  The  reaaon 
for  thla  choice  will  become  clear  later  In  the  dlacuaalon.  The  outward 
integration  from  Xjj  to  x^  la  then  carried  out.  At  x^  we  have  calculated 
a ’Teft"  value.  » "right"  value  log  derivative 


II  . 
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function.  These  two  quantities  are  equal  only  when  the  energy  parameter 
£ is  equal  to  an  eigenvalue.  The  difference 


is  a function  of  £.  (The  fector  c is  a constaid  that  will  be  specified  later. ) 
Within  some  energy  interval  around  an  eigenvalue.  D(£)  is  a well  behaved 
function  of  £ that  is  aero  at  die  eigenvalue  and  has  a positive  slope.  An 
ana-’ytically  derived  example  of  D(£)  will  be  given  later.  For  now  it  is 
suMclent  to  say  that  D(£)  is  qualitatively  similar  in  shape  to  a tangent 
function  with  roots  at  the  eigenvalues  and  a pole  between  adjacent  roots.  It 
is  somewhat  analogous  to  the  function  F(£)  defined  by  Cooley.^  The  poles 
of  D(£)  serve  as  boundary  points  that  divide  the  energy  axis  into  intervals. 
Within  each  Interval.  D(£)  is  a monotonically  increasing  function  with  a 
single  root.  As  a matter  of  convenient  nomenclature,  we  refer  to  the  energy 
interval  which  contains  the  n'th  eigenvalue  as  the  n'th  interval  and  to  die  por- 
tion of  the  function  D(£)  contained  in  this  interval  as  the  n'di  segment  of  D(£), 

A specific  eigenvalue  Is  calculated  in  a two  part  procedure.  The 
first  part  of  the  procedure  is  the  bisection  method  analogous  to  that 
described  previously.  The  Inward  and  outward  integrations  are  calculated 
and  the  poles  of  y|(x)  are  counted  [y^(x)  has  no  poles].  If  tlM  number  of 
poles  is  greater  than  n.  die  trial  energy  is  too  high, and. if  leas  than  n,lt  is 
too  low.  If  the  pole  count  Is  equal  to  n we  calculate  D(£).  If  D(£)  Is  nega- 
tive, the  trial  energy  is  below  the  eigenvalue,  and,  if  positive  it  is  above.  If 
the  trial  energy  is  low,  it  replaces  £^, and,  if  it  Is  high.  It  replaces  £^;  then 
a new  trial  energy  is  calculated' usl«v  Eq.  (6).  This  procedure  is  Iterated 
un^U  die  node  counts  for  both  £^  and  £^  are  equal  to  n, which  Indicates 
that  both  these  energies  are  wtthih  die  h'th  edrergy  Interval.  We  then 
switch  over  to  a fast  converging  numerical  method  to  quickly  locate  the 
root  of  die  n'th  segment  of  D(£).  We  have  used  the  secant  mediod^^  for 
dlls;  it  is  easy  to  program  and  has  a quite  adequate  convergence  rate.  As 
a precaution,  die  convergence  of  the  secant  method  is  monitored, and,  if  it 


■tart*  to  diverge,  the  program  returns  to  the  bisection  procedure  for  two 
more  iterations  before  trying  the  secant  method  again. 

In  order  to  obtain  a better  understanding  of  this  log  derivative 
formulation,  it  is  useful  to  examine  a simple  analytically  solvable  model 
problem.  The  simplest  bound-state  problem  is  the  infinite  square  well. 

The  potential  is  zero  in  the  interval  from  0 ^ x ^ a and  infinite  elsewhere. 

The  inward  and  outward  solutions  of  Eq.  (8),  which  obey  the  proper 
boundary  conditions^  are 

y^(x)  = kcot[k(x-a)]  (11) 


and 


y.(x)  = kcot[kx] 


The  matching  point  is  defined  by  Eq.  (9)  where 


we  can  use  the  equality  in 
this  expression  since  an  analytic  solution  is  not  restricted  to  grid  points. 
Solve  for  this  point,  substitute  into  Eqs.  (11)  and  (12), and  then  use  Eq.  (10) 
to  obtain 


I>(£)  = cktan(ka)< 


The  roots  of  this  function  occur  at  die  energy  values 


E_  * (n^^/2ua^)ti^  n*i,2,S 


which  are  obviously  the  eigenvalues  of  the  infinite  square  well  problem 


The  function  D(£)  haa  all  the  propertlet  deacHbOd  earlier.  The 


roots  of  D(E)  are  fairly  well  centered  within  their  respectlre  enerfy 
Intervals.  This  would  not  be  the  case  If  die  matching  point  for  the  Inward 
and  outward  Integrations  had  been  chosen  differently.  Let  us  assume  that 
this  point  has  some  arbitrary  value  Then  the  hinctlon  D(E)  would  be 

0(E)  > cksln[ka]/{sln[k(x^-a)]sln[kx^]}.  (16) 

It  Is  clear  that  the  roots  of  this  function  are  still  given  by  Eq.  (l$):how. 
ever,  the  poles,  which  define  the  boundaries  of  the  energy  Intervals,  are 
changed  depending  on  the  value  of  x^.  A bad  situation  occurs  If  an  Inters 
val  boundary  Is  very  close  to  a root.  This  will  Increase  the  number  of 
bisection  Iterations  required  to  establish  suitable  upper  and  lower  bounds 
before  the  secant  method  can  be  uced.  In  the  extreme  case  where  the  root 
and  Interval  boundary  coincide,  the  method  would  fall  completely. 

This  reformulation  In  terms  of  the  log  derivative  would  just  be  an 

academic  exercise  If  an  efAclent  numerical  mediod  for  solving  Eq.  (8)  did 

'*‘15 

not  exist.  Fortunately,  such  an  algorithm  Is  available.  For  convenience, 
we  briefly  restate  It  here  and  also  derive  the  formula  for  Integrating  In  the 
negative  direction  (Inward  Integration).  The  outward  Integration  formula 

Is 
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Itfam.  It  le  evident  from  the  recurrence  relation  UO)  ibat  d^  l#%  r. 

foirly  emootti  ftmctlon  of  ttie  index  n until  a critical  point  la  reached  where 
s < -1*  At  Um  next  grid  point  a^  hae  a poaltive  vdhie.  Thia  ie  the 
diecontinuity  of  the  log  derivative  function  which  occurs  at  the  pole. 
Therefore,  poles  can  be  counted  by  ^ ctaH**  procedure  of  counting. tits 
number  of  times  the  situation  < -1  is  encountered. 

Thfi  "bwmrd  inl^gm^on"  rucurslon  fo|m^  can  be  ohutned.by 
llie  sign  of  h and  also  ^hpaging  n*l  to  wf  1 In  <17)«  . ^rom 

this. we  obtain  As  formula  lor  « hy^ 
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The  Initial  term  for  fliis  inward  int^ration  is  . It  is  relateil  to  Uw  log 
derivative,  '’T  **»•  formula  i,  . 

•n  • (h/3)0{x^)l. 

We  previously  specified  this  boundary  value  to  be  y,(xjj)  • -*•  Thus, 
assuming  0(a^)  !■  finite,  a^  **  Just  as  in  the  outward  integratim^ 
formula,  total  number  of  grid  points*  musl  be  an  odd  tatogsr.  Since  ^ 

U an  even  -^ued  integer,  it^foll^  that  if  Is  also 'an  even  foieges.' 
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The  result  of  doinf  the  toward  aad  outward  iatogratloae  to  a "right"  value 
*r^*M^  and  a "left"  value  a|ia|^Mrom  which  we  calculato  ^ difference 
fuBctloD  0(S)  [aee  Eq.  <10)1  . 


D(E)  = - ■|<Xm)  - 2(h‘/3)Q(ji^) 


If  the  end  polnto  and  ^ are  located  far  ehough  laalde  toe  clasal- 
cally  forbidden  regions,  the  Initial  vaduea  sQ  and  ^e  Immaterial  sIms 
any  error  rapidly  decays  to  Insignificance  as  the  Integration  progresses 
through  the  forbidden  regions.  Under  certain  conditions,  however.  It  to 
desirable  to  reduce  toe  Initial  error  as.  much  as  possible.  An  example  of 
such  a situation  Is  when  the  elgaavalue  to  very  cldse  to  the  top  of  the  poten* 
tlal  well.  In  tols  case  tiNr  erra^  Avdayt  stowly  dsfwe  Iphigrate  torough  toe 
classically  forbidden  rsflop. . Eynstodra  gpod  hilitol  estimate,  the  Interval 
of  integrallon  required  for  toalirrer  to  become  Insignificant  to  reduced, 
tons,  savlag  computer  tone*  One  method  of  calculating  good  estimatos  of 


th«  Initial  VaTm*  Vi"  tay  tte  W,  K,  4|>iproitlin»tl«i.  Th«  W.K^  B.  *oHtio» 
Is  obtalnsd  by  asglsctlng  the  derivative  term  In  il<|»  (8)»  The  resulting 
W.  K.  B.  tlnltlel  vnlnes  ere 


The  epproprlete  values  for  Bq  and  Sj^  are  obtained  by  substituting  dkeee 
values  Into  Eqs.  (22)  and  (24).  respectively. 


The  Mimerov  mediod  Is  an  efficient  algorithm  diat  can  be  used  to 
obtain  numevleal  solutions  of  Eq.  (1^.  The  basic  formula  la  the  three 
term  recurrence  relation 


E- 


f I 


! I 


i ; 

t i 


H«r*  h la  th«  apacUif  b««waan  th«  f rl4  polnta  and  Q<xi  la  defined  by 

®q(2). 


The  renomaliaed  Nomerov  algoritbm  la  derived  from  Eq*  <50)  by 
making  two  tranalbmiationa.  The  flrat  tranafoxxnatlon^^  conalata  of 


defining  the  quantity 


*■»■(*-  \K- 


and  aubatltutlng  Into  £q.  (30)*  Thla  glvea 


F . , - U r + F , 
n+1  n n n>l 


where 


U 


k-li 


(1  - T ) *(2  ♦ lOT^). 
n n 


■ - *1  niiv  ■■ 

The  advantage  of  Ale  tranaformatlon  la  that  one  leaa  multlyllmtlon  per 


atep  la  neceaaary  to  calculate  F^  than  The  next  tranaformatlon  la 


connnonly  uaed  to  replace  a three  point  recurrence  relation  wlA  a two 
21 

point  relatloh.  Define  the  ratio 


(33) 


(34) 


(35) 


U 


■ .3 
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Subatitute  Ala  Ato  Eq.  (34)  to  obtid*  the  tera  term  recurrence  reAtion 


can  alao  be  ei^eaaed  aa  a aonthmdd  ^fraction 


Th.  la  vmIqI  aa  wa  aoon  Oum,  Wa  ipaa^r  fo  caicakte  It 

laataad  ott  bacwaa,u^UM  ttiia  qoM^tjr.  It  ^aa  not  grow  axpoMtetlaUy 
in  ttia  claaalcaily  loablddaa  ragloaa*  . 

''«*a^  ioaca  t^  Mtlal  valaa  la  ap^Ulad* 
Tha  valiMj!^ .:codltlap»a  ^(Xq)  * 0 and  4 ^ la 
R0  » •.  . Jgdta  ltfTatV  >^  Mtaay 

*aa  lipaadiataJr  asluiWi  HA 

thsaa  quaBtltl9%,aad  aalag  £ia.  (3^)  aadJW),wa  eaa  calcalala.  , 

and  X , to  wlAla  a aormaltoalloa  factor 


tonctUm  ^ ^toto  X^j. 


In  Ala  aaaaa,  wa  ara  calculating  a ranormallsad  aravafunctUm  at  aaA 
Thla  la  Aa  origin  of  Aa  daolgnation  "ranomalUad  NumaroT*'  snaAod. 

Blatt^  haa  darlvad  a laduala  ia*  coAfiatlag  |iia  darlaatlaa  of  Aa 
wavafunctlon  at  Aa  point  x^* 

aj£3g*^lft!5n  T>3(;-1»  i'^";iw.i  .;)i:oSr.  ' "■<:  . ■'•  ''‘'rii  /'.iif!  ■•'.vii  ;'.’V  f.J- ••>.(.- -i . 

~ut-arxoi  v-uTs-Ht^  'a- 

-"jjJtaq  ;‘l  .tij!.  xi'x.' 'Xtr  v. j!  rv--. 

of  ef<*r  k*  vhlek  la  poeaar  ttaa  lha  ofda^1S*'iidiTiif  St 

UnmU.  Bacaaaa  of  Ala*  Blatt^  alao  darlvad  anoAar  all^y  nora 


T-BifJ 


coinplicated  Ibxsniili^ilpr  In  trtilch  •rrof  Unn  la  of  order 

K • Sloan*  howava^*  liliii  klidwnttiHat  tlia'aiiittutiutiva  aifror  of  IIm  Nninarov 
matiiod  at  a ftxad  value  of  « la  of  order  h^.  Therefore,  It  appeara 
Xq.  (42)  oho^  probably  slve  aU  |he  accuracy  poaalbla  and  noSili^  aroold 
^ ^orm^a.  We  teated  b^  maSMte  nam* 

erlca%  on  a alinple  problem  mnd  found*  Ala  to  be  ^%aae.  Initially,  lor 
Sm  flrat  few  Iteratlona,  the  9*d»  order  medii^  gave  much  more  aoenrale 
reaulte  than  iSq;  (42).  HowaV^t  iftai^  Ihtegrafilig  «4ei^  only  a ataall 
Interval;  die  bdiddlattVe  error  Id  the’  NuAierdv  ohldtfliliatf  WaveSInctiien 
buUt  (^iad  frdtki^Sien  bd  die  alee  of  the  drroii  ld  SMa^deHeatlvea  «ma 
about  aibi^fbr  both  fbrthuldb.  Therel6re^^fe  edhcfl  efmplar 

formiiie,  Cq^  ■(42>;  f ivei  all  die  acouracy  ibaandd  MWwe  had  It; 

ihe  derivative  of  the  wavefunctloa  at  the  point  can  tie  calculated 
to  within  the  normallaatton  factor  N by  aubatltutlng  Eqa.  (39)  and  (41) 

Xq.  (42).  Then  divide  dile  1^  die  wavefup«etlon  given  by  Kq.  (40)  to 
obtain  the  log  derivative  at  [eae  »i.  (7)].  The  normallaatlon  factor  N 
cancela  In  the  numerator,  and  denominator  and  the  reault  la 


»><s>  - V 


(43) 


■fffi 


' >0 


A*  <> 


• . ■ •**■{  'A 

Although  we  are  only  conalderlng  alngle  channel  bound  atate  probtoma 
In  thla  report.  It  la  pror^  meptlonl^  reponisaHaed  Npnerov  formu< 

laa  are  eaally  generaltaed  to  a multlchanMi  matrln  formulation.  In  partle- 

??•  f m^rlx  irom  ahlch  the 

acatte^  mat^ca^  be  - 


(44) 


U'-ni  ■ :i:  J. ; 1 i 


ife  ‘ .si*;:*  ); 


'* ; i VXv 


£i 


Cqo*||^f  in  r*^rM  dlrnctlon  (U  •. « 4«er«M< 

lug  IndM  n)  »•  w«U  »•  ^ foi^s^  ^Irvctloh.  For  Um  rovorso  Itorntlmi  U 
la  conv^Bloi^  to  rodofino  tiM  ratio  [•••  Bq.  (36)) 


whoi^  lli«r<iid*  to  oaod  to  dUtliiigula)i  tia«  ratio  from  Om  ratio  {#q.  436^ 
aaad  1^'ttta  fortrord  Uorattobi  Sttbititotlav  tkt^liilo  Sq^.  (Sd)',  nr«  oUf|ii 
tiio  taro  tarm  rocorrance  formula 


Thi#lontt(da  ead  bo  lidratod  laltUd  valuo  la  t^iov^odi.^^  3^ 

valuo  iKat^orroapmiia  to  #(0^)  » 0 and  # 3 la  inn  dOrl»<^» 

vatloflNd  fommla  lortbo  tof  doriraitlvo  It  almllar  to  tba 

tloa'’0#’Bq^'-(43>.  Tlio  raaoltia-  • ^ . 


whora  la  daflnad  by  Eq.  (44).  Aa  the  Inward  Uoratlon  of  Eq.  (46)  pro- 
caoda.  tlio  raluo  of  tha  ratio  la  monltorad.  Wbaa  dia  coodttloa  a 1 
flrat  occura*  tha  ItaratUm  la  atoppad  and  dkla  point  la*  by  daflnltlon*  dM 
point  X|^.  Thla  lo  ofiproalmjataly  ^ poattl^  wbara  dM  flrat 
darlrallira  of  tba  wavalonctlon  la  aaro*  "wblch  la  dia  aama  condition  taapllad 
by  Eq.  (9).,,  Tba  outarard  Itaratlon  from  to  n^  la  carrlad  ont  nnst.  Dnr- 

locatod  botwoon  a^  and  i^|  will  giro  rlaa  to  dia  condition  < 3.  llraa* 

Ik  order  to  count  dn  nnaabar  o(^ dkicb  H^dka  Intaavnl  from  Sp  to 
atMnllor  ttmn  1 a a a M-l  and  count  dn  nombar  of  timaa  U la  nbgatlra. 

Tha  t#  Kiiiplinii  ti  aiHdtmmdrE<b.^*l63>  M*i4i|  If 

Widifb  IQi-'i  .^3  i'3,  o?  ytldsiiia  ■?;  ijiftJ 


Ths  procttdor*  outlined  abov*  liill  ba  refarrad  to  aa  mathod  A.  A 

alight  variation  of  thla  proeadara,  to  ba  rafarrad  to  aa  mathod  B,  will  ba 

daacrlbad  aaxt.  Tha  only  dlffarantai  batman  fha  taro  mathoda  la  dia  aray 

In  arhlch  dia  function  D(E)  la  calculated.  From  tha  daflnltlona  of  R and 
^ »» 

flvaarby  Sgh.  4961  aad  (dSh  rarapactkvaly,  ^ la  obvM>«*-duit  tha/Miviallty 

la  tma  only.  U E l»  an  alganvalua*  4,  foUoma  thaj^  dia  functUui^^ 


la  aaro  at  the  alganvaluaa.  With  a little  diought,  one  can  alao  ba  convinced 
diat  tha  alopa  In  ponltlva  at  the  root*  Although  thla  Juncjdon  la  pot  agunl^lo 
die  on#  datflnad  by  Sq. . It,  la  very  ahnllar  in  ahapa  and  haa  all  tha  pfOTt 
partlaa  raqulaad  for  die  oalaudatlon  of  tha  alganvaluaa.  In.  fact*  lit  mn  ha- 
ahown  that  die  taro  functlona  differ  from  each  odMr  only  by  a tairnioof  nfdaa 
h^.  Thla  foUowa  from  Eqa.  (32),  (33),  (36),  and  (45),from  which, wa  obtain 


SttbatUirttag^l^oa  Into  Bq.  (49|  aa4  in|lB||.  Eq*  ;j(^)  w|,l|h  c a h ||lyaa' 


Midtadf.A  andEpaana  fofWoalanqppily'.analW  l#Bttaid'iE.haa. ttja. aj^wnan > 

that  It  la  allghtiy  aaalar  to  program  tha  calculation  of  Eq.  (49)  ttan  Eg*.  (49Jt; 


Jbi  this  ••etUm  wm  4m»erihm  th»  r«Bonnalls«d  Nam«roir  mmtbod  for 
compotiog  •IgMiractloaa  atid  alao  darlva  a uaaliil  latarpolatloa  foranla. 

Xt  la  aaaumad  iliat  tha  alganvalua  haa  alraady  ba^jp  calcaXatajj,  Wiltk 
tha  aaargy  paramatar  aat  aqaal  to  Ihla  algaavalua,  Sq.  (37)  la  aolvad  for 
^ valaaa  of  Rj,  *2*  • • • ®9» , (46)  la  aolvad  tor  Sjj.  j.  ^^.2*  * * 

(Thaaa  qoaatltlaa  do  not  aaad  to  ba  racompatad  If  Ihay  ara  aavad  frmn  Ilia 
algaovaloa  calcolatloo. ) Tha  naat  atop  la  to  calculata  F^.  Slaea  tha  aor- 
mallaatloa  of  lito  wu'vbtoactloa  la  arbltinty,  iQit  1^  valaaa  of  iT^ 

■ for  a < M ara  caleolatad  itaratlvaly  ualag  Eq.  (36),aod  dia  valaaa  for 
' h>  M ara  calottdtad  aalag  Bq.  (43).  Ilia  Iraida  of  th#  vrdvafdactloa -at  aach 
polirt  caa  thaa  ba  compotad  aalag  Xq.  (33).  Oaa  advaalaga  tola  procadara 
' haa  ovar  a dlirdel  totogratloa  ci  toa  Nomorda  tormola^  Bq.  (30).  la^  ^t  wa 
avoid  aay  compatar  ovarfloar  problama  that  could  arlaa  daa  to  tha  aapo- 
aaatlatty.lhcraaalag  aatt^a  cd  toa  aolotloo  .la.  toa  c^a^ca4y  toii;bl4dah 
ragloaa. 

For  aoma  problama  It  la  aacaaaary  to  kaoar  tha  valua  of  toa  wava- 
. |aoctloa<al^|taU|ka  otoar  toaa  toa  avaaly  aphead  g^^ld  pcdata.  Fotr  aaanqpta^; 

. If  ora  oraat  to  latagrata  aoma  aapraaalcm  coatalalag  toa  oravafdactloo.  aoch 
aa  a matrix  alamaat  or  Fraack  Coadoo  factor,  by  Gaoaa  quadratora.  toa 

llJtii  ..Kl  - t:  ■ ’ U : ■ 0~  ^ t.  : i 'i  ' ' > f ! ■'.'S';  ■ ^ f / - i ' 

oravoloactloa  moat  ba  kaooro  at  tha  appropriate  Gaoaa  potota.  Xor  topaa 
ptdhiama  It  la  aaafol  to  hava  a aimpla  latarpolatloo  formoia  toat't^  aliW 
too  calcuiatloa  of  tha  oravafdactloo  at  ao  arbitrary  polot  to  too  aamo  accar* 
(Oqtoy  aa  It  la  ldtoBBtol^'toa  grid 'j^idta.  - 8ac|i»dl  iMhuda  caa-b*  d^lddd'Dy 
aoaaoalBg  almllar  to  toat  aaad  to  darlva  tha  Nomarov  tormola. 

Aaaopaa  that  ora  oraat  to  compata  toa^^otodo  of  lla  witoltlliiiltoi'ii  a 
poll  orhioh  la  boterooa  tha  taro  grid  polata  x,  , aad  x. 


Hi  Hliurmai  f.ai  ia.'ajV^tfCT 


it 
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trhar*  io,.  ; 

h s X,  - X,  , M»d  0 < o X 1* 

For  cettvcnlrae*, 'didllde  ''  - ^v.  - 

■T%x.  •-  •'  • ■■ 

Moko  X Txylor  oi^aiirloii  of  tbo  wavefonction  around  the  point  3^; 

♦l-l  • - «l^(4  ♦ 


(55) 


(56) 


(57) 

(58) 


ai  -,:•  '■* 

■.«*»{ 


Here  0(h^  rep^eeenCa  a term  of  order  h^;  Mo^ttply  Eq.  (5f)  liy  0**  and 
Eq.  (58)  by  a*  ^ and  add 

0-^f^  * (b0)^‘0^  + m^nw;;  + (0-b)(h^/6)ij“  + (59) 

r-  ,' i.v-  • ■ .-V.:  . ■ ' i-i 

Differentiate  Eqs.  (57)  and  (58)  twice,  add  the  reeultiaf  equatlona  and  »a^- 
tl^y  by  h^/6  to  obtain  . > , 

•:  - ■ :.  •■■■^  v - '.  :B|.  ■ :■  M : ■ , -,4  ’ -'fr* 


Subtract  4iUa>^QMIKlBi.>  (59)  la  obtain 

; i rn  fy!  i A h-u'** 


tl‘**,  +*■‘♦,.,3  - »'/*>(*!■  ♦ ♦!L,3  - M)'\ * 0>‘/*W|,  ♦ <>«»’)♦*  . <*« 

' ' ' ■ c t 

v ||Mi  oca  the  Schroadtager  aquatlaia  to  ^tmln^SM  aacaod  darlvaHira  tataaa. 

Tba  final  aj^raaelab  far  Iba  Intairpolatlan-foianula  la  ^ 
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IV,  exampijE  calculations  and  discussion 

Computer  programs  employing  all  the  meSiods  described  In  the  pre- 
ceding sections  have  been  written  and  used  to  solve  a variety  of  problems. 

Some  results  of  three  test  cases  are  presented  and  discussed  In  this  sec- 
tion. AU  resulU  were  obUlned  using  the  Aerospace  CDC  7600  computer 
and  Fortran  codes. 

The  first  two  example  calculations  are  the  Morse  potential  problem 
and  the  unsymmetrlcal  double  minimum  potential  problem  that  were  dis- 
cussed by  Wlcke  and  Harris.  ^ The  functional  form  for  both  potentials  Is 
given  by 

V(x)  » D{1  - exp[-B(x  - xj]}^  + Ae3q)[-C(x  - x^)*].  (64) 

where  the  parameter  A Is  sero  for  the  Morse  problem.  The  potential 
parameters  and  reduced  mass  are  given  In  Table  1 of  Ref.  (6)  and  also 
here  In  Table  I.  The  integration  limits  for  all  calculations  on  bo^  prob- 
Ums  were  x^^^  - l.oA  and  x^^  » 2.6A. 

Results  for  the  Morse  problem  are  given  In  Table  n.  They  are  pre- 
sented as  the  difference  between  the  numerically  calculated  eigenvalues 

26 

and  the  analytic  eigenvalues  calculated  from  the  formula 

E > l(K>0(n  1/2)  - * 1/2)^.  (65) 

n 

The  roBonnaitsed  Nomerov^^  (RN)  calculation  iM^as  done  with  2^49  gfld 
points  and  required  about  0. 73  seconds.  The  log  derivative  (LD)  method 
requires  fewer  computer  operations  per  grid  point  than  the  RN,but 
requires  more  points  to  achieve  the  same  accuracy.  These  two  effects 
almost  compensate,  m^lng  tihe  total  conqmtatlon  time  about  Uie  same. 

The  LD  results  In  Table  n were  calculated  with  2617  polntsi  the  time 
required  was  0. 79  seconds.  Calculations  to  Uils  degree  of  accuracy  can 
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Derlvatlv* 


.analytic 
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h l«  dM  fxld  spacUBf.  This  ioxtnala  Is  ossd  to  extoapplste  d»s  slfsn- 
vftluss  cslculstod  using  ssvsxsl  dUisrsnt  yuluss  of  h to  h « 0.  Results  of 
extrupoldtloa  of  RN  cslculstlons  using  261  sad  521  grid  points  sre  given  la 
Table  n.  The  total  time  required  for  dais  calculation  was  0. 2 seconds. 
Since  come  Integration  methods  discussed  In  the  literature  require  good 


Initial  estimates  of  the  eigenvalues,  It  should  he  emphaslsO^  ^t  no 


1.3.6 

► 

such  Initial  estimates  were  required  or  used  Id  these  calculation. 

A portion  of  the  fuacdOn  D(E),  calculated  lor  die  Morse  prohlsm  hy 
reaosmallsed  Numerov  method  A,  Is  shown  la  Fig.  l.  All  of  th*  » • * 


adjacent  n ■ 1 aad  a ■ 3 segments  are  plotted. 

I Q2A  Mm  and  E * 2806  cm'^  are  the  lower  aad 


Results  for  the  unsymmetrlc  double  mlalmOm  potential  proUem  are 
given  In  Thhle  III.  This  calculation  was  done  hjr  ^ renormallaof  Mumerov 
xnediod  using  2049  points.  The  time  required  was  0. 86  seconds.  These 
values  differ  by  a small,  but  significant,  amount  from  those  oUalned  by 
Wlcke  and  Harris  [see  Table  H la  Ref.  (6)].  BehSuse  of  this  dlHsrhne  we 
checked  opr  results  by  doing  sn  Independent  calculation  with  a bqsls  ^no- 
tion eagiMif  Ion  method  udllaUgi  a 60  term  harmonic  oscillator  fimctlon 
eaqianslon.  The  agreement  was  encelleat,  convlaclag  us  daat  das  results  la 
Table  IB  ere  correct.  The  dttfereace  between  bur  results  and  dlOihdf 
Wlcke  a^  Harris  Is  prohahlf.dpe  to  differences  la  the  calcuUteg^^a^lpes 
of  Wm— ^ B and  i«a*d  hi  Am  respective  cals. . 

eulatlcas.  la  order  •dhaMMdteb^aoas  adWhtdlss  dtodM^^l^^ 

solvlag,  ere  give  our'ftila'S’af-'ll  4a  lddedi*ei^dlgltif^v-4l'«'  l**8dtBWI^i^^8ji8| 
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he  dowe  meee  econoialoaUy  by  ualaCithf  Blchairdson  ea^apolptloia  tech- 
nique. **^'^*r  It  h^  been  vSrUled  by  nuDasprlcal  calc«tlattonf.that  dw  trunr 
cation  error  for  both  the  RN  aad  LD  methods  Is  given  by 


Table  m.  Elfanealuea  lor  ti»a  Unaynanatrlc  and  Symmatrlc 
Double  Minimum  Potential  ProMema 
Calculated  ualng  Uie  Renormallaed  Numerov  MedMd  A. 
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0 

1202. 500 

i 

0.128fll928 

1 1 

3205. 207 

0.128811949 

i 

2 

4227.329 

0.409028941 

2 

9144.251 

0.405920240 

4 

6064.241 

0. 650944055 

5 

7092. 679 

^651100997 

1 

6 

7614. 622 

0.864917277 

7 

8911.545 

0. 872946249 

l' 

8 

9095.696 

1.01722996 

: 

. j 

9 

10208.290 

1.07895209 

i 

1 
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10 

' 10869.289 

1. 18917992 

( 

1 

i 

11  u 

11482.479 

1.20119279 

12953.^ 
t2fn.479  ' 
lH90.4iS 
14429. 9S0 
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1.42924820 
1. 99718929 
1.9H40809 
1.84277820 
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TImt  foBctioa  D(EX  «««  -eomputod  for  thU  prpbl«m.«oi  U WM  ob«f  viw4 
that  tha  ahivd  aaargir  Uitetval  la  vai^  Aanraar*  TImi  lowar  hatiadarr  ^*^.  . 

locatad.a*E  ••  5097  and  tha  appar  haaadhMy  la  aiE  • ildl  caa^f. 
Bacaaaa  thla  Interval  la  ao  narrow*  more  blaactloa  Itaratlona  Otaa  naual 
are  reqnlrad  before  we  awltch  over  |o  the  faat  converging  aeoant  SM^bod, 
The  email  wldtti  la  due  to  ^e  condition  that  ^ ntatchlng  joint  la  located 

near  the  outer  minimum  w**ile  the  n « 3 eigenfunction  la  mainly  localliaed 
around  the  Inner  minimum.  In  order  to  demonatrate  thla*  we  altered  ttw 
program  ao  that  the  matching  point  waa  located  near  die  Inner  inlnlmnm. 

In  dila  caae  the  third  energy  Interval  waa  much  wider  with  die  lower 
boundary  at  E « 4230  cm~^  and  the  vqiper  boundary  at  E e 6030  cm*^*  A 
more  aophlatlcated  program  could  be  written, which  would  automatically 
locate  the  matching  point  In  tha  region  where  the  wavefunctlon  la  hutlnly 
locallaed}  however*  we  did  not  do  thla*  It  ahould  ba  emphaalaed  diat  ho 
matter  how  narrow  an  Interval  la*  our  procedure  never  falla  to  converge 
to  the  dealred  eigenvalue;  It  ahnply  tequlrea  a few  more  blaectlon  ttera> 
tlona.  Thla  la  die  reaaon  that  die  double  minimum  pmhlem  req>nlred 
allghtly  more  time  to  aolve  than  die  Morae  problem.  In  a more  extreme 
teat  caae  we  changed  the  Qauaalan  amplitude  factor  In  the  po^tlal  to 
A m 20000*  In  dila  caae  the  n * 3 energy  Interval  waa  ao  extremely  i^irrow 
that  convergence  waa  achieved  entirely  by  the  blaectlon  medio^  the  ae<mnt 
method  waa  never  called  Into  uae* 

The  third  example  calculatloa  la  a aymmetrlc  double  mlnlmnm  peten- 
tlal  problem*  The  potential*  In  arbitrary  unite*  la 


Numarov  medMd  A*  are  given  In  Table  HI*  Th^^  ht^gre^^JU^a^  ^ ^ . 
42.0*  the  number  of  grid  polnta  waa  4001*  and  the  time  re^nlifed  waa  1*9 
aeceada.  Maee  the  lower  elgenvalnea  are  donhleta  wldi  very  mull 


•plKttikgSf  eonv«fv«ttft4i^lo  OuMM  valttM  will  b*  «iaaoat  AatlMly  bf  tiw 
bl««ctlott  pco«44ttf*.  This  fbet  In woBjimaCira  wltb  tba  tevf*  MAbar  bl 
tntbf  ratUm  pbbUrt  accowMw  for  tiMr  lacraaaad  tlma  lo  adva  ftta  paobtanu 

Anotiltar  mallwd'^  a^lcb  takaa  account  ot  tbe  ayinmatry,  can  ba  uiia4 
to  aolva  ^ia  probiam  more  aitflclantljr*  The  tntaf  ration  la  canlad  out 
over  only  ibatf  tiia  r«n|a/  ® *toa»  * the  aymmalrlc  a^ 

antlaymmetrlc  tolutlona  are  obtained  la  aeparate  calculatlona.  Thla 
metibod  la  mere  efficient  for  two  reaaonat  The  Integration  coeera  only  half 
the  range*  and  alnce  the  elgenvaluea  no  Icnger  appear  aa  doubleta*  fewer 
Ueratlona  by  &e  blaectlon  method  are  required.  Ibe  aattaymmetrlc 
elgenfunctloaa  are  calculated  by  merely  clwngtng  ibe  atartlng  point  of  the 
outward  Integration  from  Xq  « >2  to  « 0.  Since  the  waveflianctlbn  la  aero 
at  d&la  point*  the  Inner  bou^ary  condition  la  the  aaine  aa  have  dlacuaaed 
prevlooaly.  The  aymmetrlc  wavefunctlona*  on  the  other  hand*  have  herb 
alope  at  X ■ 0.  From  Ihla  we  obtain  * 6,  vdilch  la  auhatituti^  Idb 

Bq.  <22)  to  calculate  dw  Initial  term  a^  of  the  log  derlvatl^^  lolufiba.  The 
Initial  term*  R^*  of  a renormallaed  Numerov  aolutton  la  alao  eaay  to 
derive.  Let  p^|*  00*  and  be  die  valuea  of  the  wavefunOtlon  at  die  grid 
polnta  x^i  « •k*  M 0*  and  x^  ■ h.  From  aymmetry  It  folidwa  thit 
0^1  ■ 0|  tod  alao  F^i  ■ F^.  Subatltuttalig  dila  Into  Bt|.  f24)  and  ualng 
Xq.  ffd)  we  obtain  die  dealred  Initial  term 


where  Uq  la  defined  by  Sq.  (35). 

A calculatloa  of  the  aymmetrlc  and  antllymmetrlc  aoluttoaa  waa  cal- 
culated dda  way  uatag  the  reaoamallaed  Numerov  mediod.  The  reaulta  are 
exactly Mate ’ad  dyatilift'prev^^  aad'Uated  (iii^Tayeiffl.  '^'l^  mm-*  ' 
bar  of  grid  iiiiala  niaed  daa'Sdll'  andT  llte  latal'  thm'  ireqiAred  %o  cMtolale' 
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